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SIMPLIFIED ANALYSIS OF NUCLEAR FUEL PIN SWELLING 


by Armln F. Lietzke 
Lewis Research Center 

SUMMARY 

The effect of fuel swelling on strains in the cladding of cylindrical fuel elements is 
analyzed. Simplifying assumptions are made to permit solutions for strain rates in 
terms of dimensionless parameters. The most significant assumptions included in the 
analysis are (1) the fission gases appear as bubbles dispersed throughout the fuel, (2) 
the bubble gas pressure is in equilibrium with the creep strength of the fuel, (3) uniform 
temperature exists in the fuel and in the clad, and (4) no axial force is transmitted be- 
tween the fuel and the clad. 

The results of the analysis are presented in the form of equations and graphs of di- 
mensionless parameters which illustrate the volumetric swelling of the fuel and the 
strain rate of the fuel pin clad. A parameter is defined in terms of operating conditions, 
the value of which determines the strain rate of the clad relative to the volumetric swell- 
ing rate of the fuel. 

Although the assumptions used to simplify the analysis restrict its application to 
special cases, the results show the relative importance of operating conditions and 
material creep properties for these cases and give some insight as to their importance 
for other cases. The assumptions involved make the analysis most suited to long oper- 
ating lifetimes with low burnup rates, such as may be anticipated for space power reac- 
tor systems. 


INTRODUCTION 

Fuel elements of nuclear reactors are subject to internal stresses which may lead 
to large dimensional changes or gross failure. Of particular significance for long life, 
high power density reactors are those stresses resulting from a high concentration of 
fission products. Fission products cause the fuel to swell if retained in the fuel or, if 
released from the fuel, the gaseous products exert pressure on the clad material. In 
either case, some means must be provided to accommodate the fission products, other- 


wise an intolerable stress will be produced in the clad. 

The fuel element geometry frequently considered for severe operating conditions is 
a cylindrical fuel pin incorporating a void space, as illustrated in figure 1. The cylin- 
drical clad forms an efficient pressure container, and the void space accommodates any 
fission gases released from the fuel and some degree of fuel swelling. 

Several computer programs have been written to calculate stresses and strains in 
cylindrical fuel pins (e. g. , ref. 1). Computer programs which take into account all as- 
pects of the design problem, including heat transfer, thermal stress, and property vari- 
ations as well as fission product swelling, are necessarily very complex. Design calcu- 
lations utilizing such computer programs may require extensive iteration to arrive at a 
satisfactory solution and might best be utilized to check final designs. 

The purpose of this report is to investigate only those stresses and strains caused 
by fuel swelling and through this simplification to express the calculated strains in 
terms of parameters that can more conveniently be used for design purposes. The 
swelling problem treated herein is the case where all fission products, both solid and 
gaseous, are contained within the fuel and contribute to fuel swelling. The analysis ne- 
glects elastic strains and first-stage creep. The analysis involves (1) the treatment of 
fuel swelling, and (2) the interaction of fuel and clad as the clad restrains the fuel from 
swelling radially outward. 

Because fuel swelling is not fully understood, a calculational model must be hypo- 
thesized to explain the mechanism of the swelling process. Several models have been 
proposed. A model was proposed by Wyatt (ref. 2) to calculate swelling in uranium 
metal. This model assumes the gaseous fission products to accumulate at nucleation 
sites to form bubbles dispersed throughout the fuel material. The gas pressure in these 
bubbles stresses the fuel material, which strains in accordance with its creep proper- 
ties. Essentially the same model was later used by Foreman (ref. 3) to calculate swell- 
ing rates in uranium. Foreman, however, used a more rigorous treatment of the stress 
distribution in the fuel. The methods of reference 3 are followed in the present work to 
calculate the volumetric swelling of the fuel. 

A model for fuel swelling was also proposed by Barnes (ref. 4), which assumes the 
gas pressure forces to be in equilibrium with the forces of surface tension. The creep 
strength of the material was neglected. The model of reference 3 is preferred for the 
work herein for two reasons. First, high temperature applications will probably employ 
refractory fuel compounds having high creep strength. Secondly, although surface en- 
ergy effects are important initially, long operating times at high temperature will prob- 
ably result in bubbles large enough that surface tension effects may be neglected during 
most of the operating time. 

The interaction of fuel and clad is analyzed by the use of multiaxial creep equations 
applied to thick-walled cylinders. The assumption is made when applying the boundary 
conditions that no axial force is transmitted between fuel and clad. 
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The analysis results in a set of equations expressing the strains in terms of creep 
properties, operating conditions, and geometric dimensions. Graphs of some of the 
equations are presented to aid in computations. A numerical example is given to illus- 
trate the use of the figures and equations. 

The equations to perform the numerical integrations included in this work were pro- 
grammed by J. L. Thompson. 


ANALYSIS 

The methods used to calculate stresses and strains in a cylindrical fuel pin are as 
follows: Fuel volumetric swelling rates were calculated according to the model proposed 
by Wyatt (ref. 2) and modified by Foreman (ref. 3). These volumetric swelling rates 
are assumed to be independent of clad restraint, that is, clad restraint affects only the 
distribution of fuel strain in the radial, tangential, and axial directions. Using these 
volumetric swelling rates, the fuel and clad were then analyzed as radially loaded thick- 
walled cylinders having a common pressure at their interface; the fuel being radially 
loaded on its outside surface and the clad being radially loaded on its inside surface. 

Multiaxial stresses appear in the fuel swelling model and in the fuel and clad cylin- 
ders as a result of their mutual interaction. The analytical treatment of steady-state 
creep under multiaxial stress appears in the prior literature (e.g. , ref. 5). The equa- 
tions based on this prior treatment are given in appendix B. All symbols are defined in 
appendix A. 


Fuel Swelling 

The analysis of fuel swelling caused by fission gas containment in the fuel follows 
that given in reference 3. The calculational model assumes gaseous fission products to 
accumulate at nucleation sites to form bubbles dispersed throughout the fuel. The fuel 
material is divided into polyhedral cells, each of which contains a gas bubble at the cen- 
ter. Each cell is approximated by a sphere of equal volume. The fuel swelling problem 
for gaseous fission products is thereby reduced to a determination of the creep of an in- 
ternally pressurized thick-walled sphere. This latter case is also treated in reference 5. 

The model neglects surface tension and assumes the fission gas pressure to be in 
equilibrium with the creep resistance of the fuel material. The only justification for 
this assumption other than to simplify the analysis would be to show experimentally that 
the fission gas bubbles are indeed large enough to make surface energy effects negligibly 
small. Until such evidence exists for any fuel in question, it may be argued that while 
surface energy effects are important initially (when the bubbles are small), the bubbles 
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will grow by creep of the fuel and/or coalescence of bubbles. Long operating time will 
probably result in bubbles large enough that surface energy effects may be neglected 
during most of the operating time. Therefore, the assumption is expected to be most 
applicable for long operating times, which is the major concern of this report. 

Following reference 3, the fractional change in fuel volume caused by gaseous fis- 
sion products x for a constant fission rate b/r is given by 

o 



Equation (C23) can be derived from the multiaxial creep equations of appendix B. This 
derivation, shown in appendix C, is included in this report for completeness. Equation 
(C23) is equivalent to the results of reference 3, except that in reference 3 x g was as- 
sumed to be small compared to 1 and the term (1 + x ) does not appear. 

& 

Swelling caused by solid fission products is assumed, as in reference 2, to be di- 
rectly proportional to the number of fissions based on atomic and molecular volumes. 
This assumption may be expressed as 


C bt 

x = - 5 - (C24) 

r 

The total change in fuel volume expressed as a fraction of the original volume is 
given by 


x = Xg + x g (C25) 

The total volumetric swelling rate of the fuel is required for analysis of the plastic 
deformations in the fuel and clad cylinders. The engineering bulk strain rate for a con- 
stant fission rate is given by 


x = Xg + x g (C26) 

where x and x may be obtained from the derivatives of equations (C23) and (C24). 
Thus, 
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(C22) 



b 

r 


(C27) 


The strain rates are more accurately expressed on the basis of true strain rather 
than engineering strain. The true bulk strain rate is related to the engineering bulk 
strain rate according to 


v = -5_ (B18) 

1 + x 

and therefore 

v = ln(l + x) (B19) 

The engineering strain rate x and the true strain rate v are essentially equal for 
small values of x. 

Combining equation (B18) with equations (C25), (C26), (C22), and (C27) yields 


'3C y KT f Nbt\ n 
nT 


v = . 


n-ll 



in 


+ C, 


1 + x„ + C, 


bt 


(C28) 


Stresses and Strains in Fuel and Clad Cylinders 

Equations for stress and strain in the fuel and clad cylinders are derived in appen- 
dix D. The derivation involves the following assumptions: 
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(1) The fuel and clad operate at uniform but possibly different temperatures; that is, 
radial temperature variations within the fuel and clad are neglected. 

(2) The cylinders are long compared to their wall thickness, so that plane strain 
exists in the axial direction; that is, axial strain is independent of radius. 

(3) The hydrostatic pressures at the inside surface of the fuel cylinder and the out- 
side surface of the clad cylinder are zero. 

(4) No axial force is transmitted between the fuel and clad cylinders. 

Applying these assumptions, the tangential strain rate at the fuel-clad interface is 
given by the following equation from appendix D: 


where 




1 = 0 


a = 


n f/ n r 


3\ 1- (nf/ n c) 



(D30) 


(D31) 


Therefore, the parameter a determines the relative tangential strain rate at the fuel- 
clad interface (3c tm/v)- 

The quantity j3 is a function only of the radius ratio <p a or cp-^ for a given value 
of n. Values of j3 have been calculated by numerical integration (see appendix D) and 
are shown plotted in figure 2 for several values of n. 

The relative tangential strain rate at the fuel-clad interface (3e tm /v) can be obtained 
from equations (D30) and (D31) and figure 2. All other strain rates can be determined 
from the following equations written in terms of the strain rate at the interface: 



(D33) 
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(D34) 
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r 


= I + 
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(D35) 


When applying these latter equations to the clad, v is taken to be zero. Values for 
the constant A have been calculated by numerical integration (see appendix D) and are 
shown plotted in figure 3. 

The strain rates are, in general, functions of time and must be integrated with re- 
spect to t to obtain the strain in the time period from time zero to time t according to 


e = 



(D36) 


It is desirable to express the strains as engineering strains to give a measure of the 
change in dimension relative to the original dimension. True strain e can be converted 
to engineering strain 6 according to the following relation from appendix B: 


6 = e e - 1 


(B13) 


RESULTS AND DISCUSSION 

Strains in the fuel and clad of a nuclear reactor fuel pin caused by fuel swelling can 
be calculated from the equations developed in the appendixes and presented in the analy- 
sis. These equations and their application are discussed in the following sections. 


Fuel Volumetric Swelling 

The volumetric swelling of the fuel caused by gaseous fission products is given by 
equation (C23). This fractional change in volume x^ is seen to be a function of the fuel 
characteristics (N, a^, and n^), the design conditions (T, b, and r), and the operating 
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time t. The solution to equation (C23), obtained by numerical integration, is presented 
in figure 4 for several values of n. The amount of swelling increases exponentially with 
time, and depending on the fuel and operating conditions can reach appreciable values. 
Volumetric swelling of the fuel caused by gaseous products should at least be limited to 
the point where the gas bubbles become interconnected leading to fission gas release 
from the fuel. Barnes (ref. 4) estimates that this point is reached at a volumetric swell- 
ing of 33-^ percent. Fuel swelling probably ceases at this point. 

The amount of swelling given by equation (C23) and figure 4 is based on a particular 
swelling model. The model is subject to question, as would any model proposed on the 
basis of present knowledge. The fact that surface energy effects were neglected in the 
model implies that the analysis is not well suited to short-term irradiations when the 
fission gas bubbles are expected to be small. Neglecting surface energy effects is felt 
to be of minor consequence for long irradiation periods but leads to conservative designs 
for all cases. Swelling models which include surface tension must incorporate one or 
more additional assumptions not required here; namely, the mechanism of bubble coales- 
cence and/or the bubble density. As better swelling models become available, they could 
be applied to the stress-strain analysis of the fuel and clad cylinders (appendix D) which 
is independent of the swelling model. 

Swelling predictions of equation (C23) and figure 4 involve the assumption that volu- 
metric swelling is unaffected by clad restraint. Actually, clad restraint would be ex- 
pected to decrease the volumetric swelling of the fuel (see appendix C). However, the 
presence of a void space inside the fuel pin may allow the volumetric swelling of the fuel 
to approach the unrestrained swelling. In any case, the assumption used herein leads to 
a conservative design. 

The amount of swelling depends on the fuel operating temperature, as can be seen 
from equation (C23). The fuel temperature T f and the coefficient a in the creep rate 
equation (which is a function of temperature) both appear in equation (C23). The assump- 
tion of uniform fuel temperature could lead to significant errors where large temperature 
variations exist. Therefore, the analysis presented herein is not valid for such cases. 
Treatment of fuel pins where fuel swelling is a function of radius requires that the analy- 
sis be programmed for a digital computer and is beyond the scope of this report. 

The amount of swelling caused by gaseous fission products increases with tempera- 
tures and is the major contribution to fuel swelling for high temperature applications. 
However, swelling caused by solid products is not negligible and its contribution (eq. 

(C24)) is included in the total swelling given by equation (C25). 

Recent information obtained in the United Kingdom indicates that not all of the gas- 
eous products participate in bubble growth because some of the gaseous products may 
exist in a supersaturated solid solution with the fuel material under irradiation condi- 
tions (ref. 6). This unexpected effect is also under investigation as part of the fuel 
swelling studies at the Battelle Memorial Research Institute (ref. 7). A knowledge of 
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the amount of this increased solubility could be incorporated in the present analysis by 
using a properly reduced value of Cy. Some evidence also exists to indicate that the 
creep strength of the fuel material under irradiation may be lower than that measured 
out-of-pile (ref. 7). This difference in creep strength has been attributed to local over- 
heating of the fuel in the vicinity of the fission sites. 


Strains in Fuel and Clad 

The volumetric swelling rate v at any time t can be calculated from equation 
(C28) once the value of Xg at time t is known. The tangential strain rate at the fuel- 
clad interface may then be determined from equations (D30) and (D31) using figure 2 to 
obtain values for /3 in equation (D31). Equation (D30) is shown graphically in figure 5. 
The relative tangential strain rate (3e tm /v) at the interface is a function of a and 
nj/n c . At small values of a (representing a strong fuel and a weak clad), the relative 
tangential strain rate approaches a value of 1. That is, the tangential strain rate is 
nearly the same as if the fuel were unrestrained and equal to one-third the volumetric 
swelling rate. At large values of a (representing a weak fuel and a strong clad) the 
relative strain rate approaches zero. That is, the clad strain rate is nil and the fuel is 
forced to strain axially and inward into the central void space if such a space exists. In 
between these two limits, as a approaches a value of 1, is a region where the strain 
rate is more sensitive to fractional changes in the value of a. 

The value of a is dependent on the value of v which is a function of t. Therefore, 
the strain rate is also a function of t and, in general, must be integrated (either nu- 
merically or graphically) with respect to time to obtain the strain in the time period 
from time zero to time t, according to equation (D36). If, however, v does not vary 
greatly with t or if n^/n, is nearly equal to 1, then a will be insensitive to time. In 
such a case, the quantity 3c^ m /v will be nearly constant, and no time integration is 
necessary. The strain ej. m at any time t for these cases may then be evaluated with 
little error from the values of 3e tm /v and v/3 at time t as given by 



The magnitude of the error introduced by this equation obviously depends on the 
variation in the value of 3e^ m /v with time and may be evaluated from the following 
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Thus, if 3e tm /v is constant, the second term is zero. If 3e tm /v varies with time, 
the maximum error introduced by neglecting the second term is 

v 

max 
3 



Because v increases with time, the maximum value of v occurs at time t and the 
maximum error introduced by neglecting the second term is 


v 

3 



Strain rates at other locations in the fuel and clad cylinders may be calculated from 
equations (D33), (D34), and (D35) using figure 3 to obtain values for A. Again, if v 
does not vary greatly with t or if n f /n c is nearly equal to 1, no time integration is 
necessary, and the strain at any time t may be evaluated from 



A conventional form of the creep rate equation was used in the analysis. There are 
limits which should be imposed on any equation of this form for practical materials. 
Limits on strain are obviously necessary because every material will ultimately rupture. 
Most designs will have strain limits on diametrical swelling of the clad dictated by other 
considerations, such as heat transfer. Limits on stress may also be necessary to re- 
strict the range over which any one creep rate equation applies. It may therefore be ad- 
visable, if high stress levels are suspected, to check the stress level using equations of 
appendix D and choose a strain rate equation in accordance with the calculated value of 
stress. 

The equations resulting from the analysis which predict the creep rates in the fuel 
and clad are independent of fuel pin size. Creep rates are seen to be functions of terms 
which are entirely dependent on dimensionless parameters defining the operating condi- 
tions, material properties, and radius ratios. 
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Although the assumptions used herein restrict the analysis to special cases, the re- 
sults show the relative importance of operating conditions and material creep properties 
for these cases and give some insight as to their effect for other situations. The as- 
sumption of uniform temperature in the fuel and clad can never be achieved in practice 
and can only be approached with small diameter fuel pins, high thermal conductivity ma- 
terials, and low burnup rates. Reactor applications requiring long operating lifetimes, 
such as space power systems, imply low burnup rates. These are the applications of 
most interest here. Such applications coupled with a fuel of high thermal conductivity 
will operate with a small temperature drop across the fuel. Other reactor applications 
which operate with a large temperature drop across the fuel and/or large variations in 
burnup across the fuel require that the analysis be programmed for a digital computer. 


Numerical Example 

A numerical example is presented to illustrate use of the equations and graphs and 
also to illustrate typical results obtained from the analysis. The following inputs were 


selected for this example: 

Absolute fuel temperature, T^, K 1500 

Coefficient in equation (B4), for fuel, a^, ^N/cm^ ^ (hr 1.026x10 

Exponent on stress in equation (B4), for fuel, n^ 4.0 

3 22 

Total uranium atoms per unit volume of fuel material, N, atoms/cm 3.3x10 

Fractional change in atomic or molecular volume of solids occurring 

as result of fission, C„ 1.2 

b 

Absolute clad temperature, T , K 1400 

Coefficient in equation (B4), for clad, a c , ^N/cm^j c (hr *) 3.34x10 

Exponent on stress in equation (B4), for clad, n c 3.7 

Fraction of fission products in gaseous state, Cy 0. 125 

Boltzmann constant, K, N-cm/K-mol 1.381x10 ^ 

Burnup fraction, b 0.03 

Total irradiation time to achieve burnup fraction b, r, hr 20 000 

Fuel radius ratio, <p a 0.3 

Clad radius ratio, 1.2 


The volumetric swelling caused by gaseous fission products x may be obtained 

o 

from figure 4. The value of the abscissa in figure 4 is 


11 


3 Xl. 026X10” 18 x 20 000 

(3 x0. 125 X 1. 381X10" 21 x 1500 X3. 30xl0 22 x0. 03 \* 

2(4 + 1) 1 

l 4 )\ 


1/(4+ 1) 


xl = 0. 0965 — 
r r 

The values of t/ r at various values of time t are given in table I along with the values 
of Xg read from figure 4. 

The volumetric strain rate of the fuel is calculated from equation (C28). For ex- 
ample, at t = 10 000 hours, 


3x1.026x10~ 18 /3xQ. 125 Xl. 381xl0~ 21 xl500 X3. 30xl0 22 x0. 03 xlO 000 ^ 
2 V 4x20 000 i 


0. 059 


4-1 


1 - 

/ 0.059 X 1 / 4 ! 


\1 + 0.059/ 


1.2x0.03 
20 000 


V = 


1 ■ 0.059 i (1-2)(0.03)(10 000) 
20 000 


= 1.02X10 -5 


Values of j3 are obtained from figure 2 : 


/3j = 0. 794 at cp a = 0. 3 and n^ = 4 

j8 c = 0. 159 at (p^ =1.2 and n c = 3.7 

Values of a are computed from equation (D31). For example, at t = 10 000 hours, 
a 1.026x10" 18 / 3 \ l-(4/3.7) / q 159 

( 3 . 34X10- 20 ) 4/3 ' 7 \1-02X10-V V°' 794 

The relative tangential strain rate at the fuel-clad interface 3cj m /v is read from 
figure 5 or calculated from equation (D30), as given in column 6 of table I. The tangen- 
tial strain rate at the interface is 


J 1 = 0. 727 
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For example, at t = 10 000 hours, 


e tm “ °- 593 


r 1.02x10 


-5 N 


2.02X10" 6 hr -1 


Other strain rates are calculated from equations (D33), (D34), and (D35), noting 
that for the clad (1 < <p < < p-^) v is taken to be zero. For example, A c = 0. 397 (from 
fig. 3(b)), then at t = 10 000 hours, 


(0 - 2. 02x10” 6 ) (1 + 


0.397x1.2" 


e tb= 0 


Vi 


= 1.52X10” 6 hr" 1 


1 . 2 ^ 1 + 


0. 397 

Vi 


The strain rates at various positions in the fuel and clad vary with time. The true 
tangential strain rate at the fuel-clad interface is shown in figure 6 as a function 

of time. Integrating this function from time zero to time t either numerically or 
graphically yields the strain at any time t (eq. (D36)). For this numerical example, 
3kt m /v does not vary greatly with time (see table I). Therefore, no time integration is 
necessary, and the approximate equations given in the section RESULTS AND DISCUS- 
SION are applicable. Thus, at t = 10 000 hours, 


€ tm " 




+ x. 



3 


°- 5 - 9 - 3 lnfl + 0.059 + (1. 2)(0. 03)(0. 5)1 = 0.015 
3 ' 


Engineering tangential strain is then determined from equation (B13) 

6. = e °' 015 - 1 = 0.015 

tm 

Engineering tangential strain as a function of time is shown in table I and in figure 7 
for three locations: the fuel-clad interface, the clad outside radius, and the fuel inside 



radius. The axial strains 5 for the fuel and clad are also shown. All of the strains 

z 

increase monotonically with time. It can be seen that the clad increases in diameter at 
both its inside and outside surfaces, whereas the fuel increases in diameter at its out- 
side surface but decreases in diameter at its inside surface. The tangential strain in 
the clad is highest at its inside surface and equal to 4. 1 percent at the end of the as- 
sumed life. The change in outside diameter of the clad is 3. 0 percent at the end of life. 

The axial strain in the fuel is the largest of all the strains. The fuel increases in 
length by 11.4 percent. The clad decreases in length by 1. 5 percent. In order to ac- 
commodate the relative axial expansion of fuel and clad, an axial void space of at least 
12.9 percent of the fuel length must be provided. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 8, 1969, 

120-27. 
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APPENDIX A 


SYMBOLS 


A grouping of constants defined by 

eq. (D13) 

a coefficient in creep rate equation 

(eq. (B4)) 

B grouping of constants defined by 

eq. (D17) 

b burnup fraction, fraction of total 

uranium atoms fissioned in 
total time r 

Cj,C 2 constants of integration 

C g fractional change in atomic or 

molecular volume of solids oc- 
curring as result of fission 

Cy fraction of fission products in 
gaseous state 

K Boltzmann constant 

l arbitrary length 

N total uranium atoms per unit 

volume of fuel material 

n exponent on stress in creep rate 

equation (eq. (B4)) 

p pressure 

r radius 

T absolute temperature 

t time 

u rate at which an incremental 

radius dr is changing 

V fuel volume at any time t 


v true bulk strain caused by solid 

and gaseous fission products, 
defined by eq. (B17) 

v true bulk strain rate caused by 

solid and gaseous fission 
products, defined by eq. (B15) 

x fractional fuel volume change 

caused by solid and gaseous 
fission products, engineering 
bulk strain defined by eq. (B16] 

x fractional fuel volumetric swell- 

ing rate caused by solid and 
gaseous fission products, engi- 
neering bulk strain rate de- 
fined by eq. (B14) 

Xg fractional fuel volume change 

caused by gaseous fission 
products 

Xg fractional fuel volumetric swell- 

ing rate caused by gaseous 
fission products 

x_ fractional fuel volume change 

caused by solid fission 
products 

x fractional fuel volumetric swell- 

ing rate caused by solid fission 
products 

a defined by eq. (D31) 

/3 defined by eq. (D32) 

y defined by eq. (Bl) 
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Subscripts: 


engineering strain, defined by 
eq. (Bll) 

6 engineering strain rate, defined 

by eq. (B9) 

e true strain, defined by eq. (B12) 

e true strain rate, defined by 

eq. (BIO) 

e* effective multiaxial unit strain 

rate 

a stress 

a* effective multiaxial stress 

r total irradiation time to achieve 

burnup fraction b 

(p radius ratio, r/r m 


a inner radius of fuel cylinder 

b outer radius of clad cylinder 

c clad 

f fuel 

i inner radius of fuel sphere (gas 

bubble radius) 

m reference radius, radius at fuel- 

clad interface (cp = 1) 

o outer radius of fuel sphere 

r radial 

t tangential 

z axial 

0 original value at t = 0 

1, 2, 3 components in triaxial stress 

condition 
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APPENDIX B 


MULTIAXIAL CREEP EQUATIONS 

The analytical treatment of steady-state creep under multiaxial stress appears in 
the literature (e.g. , ref. 5). The relation between stress and strain rates in this prior 
treatment is based on the assumption that at any point in a stressed body the principal 
shear strain rates are proportional to the principal shear stresses 



(Bl) 


where e is the strain rate, o is the stress, y is a proportionality factor to be evalu- 
ated, and the subscripts 1, 2, and 3 refer to the three mutually perpendicular strain 
axes in a multiaxial system. 

The sum of the strain rates is equal to the rate of change in volume. Therefore, 


e l + e 2 + e 3 


(B2) 


The change in volume of a solid under plastic deformation is usually taken to be zero 
(v = 0). The more general form (eq. (B2)) is used herein so that it can be applied to the 
fuel cylinder (solid plus gas bubbles) which is increasing in volume as fission products 
are generated. 

Combining equations (Bl) and (B2) yields 


_ v + y 
3 3 

_ v + y 
3 3 

= v + y 
3 3 


(2a 1 - a 2 - cr 3 ) 

(B3a) 

( 2ct 2 " °3 " CT l) 

(B3b) 

(2a 3 - a 1 - a 2 ) 

(B3c) 


A creep rate equation of the following form is assumed for a given temperature in a 
multiaxial system 


where the effective stress and strain based on distortion energy are defined by 



The quantity y in equation (B3) is determined from a uniaxial stress condition by 
setting Og = Ug = 0. Then, eg = f 3* and equations (B2), (B5), and (B6) reduce to 





1 


v 

3 


a* = Oj 

Applying this uniaxial stress condition to equation (B3a) and substituting from equation 
(B4) yields 


r = 3a a *n-l 
2 


(B7) 


Combining equations (B3) and (B7) yields 


*n-l 


= I + i l _(2 ff l‘ a 2“ ff 3 ) 

(B 8 a) 

• __*n-l 

3 2 

(B8b) 

• _ _*n-l 

= 3 + ~2 (2<T ^ " a l ~ a 2 } 

(B8c) 


Thus, equation (B8) is the relation between stress and strain in a multiaxial system. 
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The foregoing equations are strictly applicable for true strain rates. It is desir- 
able, however, to express the strains as engineering strains to give a measure of the 
change in length or volume relative to the original values. Defining engineering strain 
rate as 


5 = 


dZ 

Z 0 dt 


and true strain rate as 


e = 


dZ 
Z dt 


then the engineering strain is given by 


/ 


5= / 


l o l o 


and the true strain is given by 


e = 


J, 1 l » 


Therefore, the engineering and true strains are related according to 


e = ln(l + 6) 


or 


<5 = e e - 1 


Similarly, the bulk strain rates are defined as 


x 


dV 

V 0 dt 


V 


dV 
V dt 


(B9) 


(BIO) 


(Bll) 


(B12) 


(B13) 


(B14) 


(B15) 
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Therefore, 


and 


V - V, 


v = In — 


1 + x 


v = ln(l + x) 


(B16) 

(B17) 

(B18) 

(B16) 
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APPENDIX C 


FUEL SWELLING 

The analysis for fuel swelling caused by fission gas containment in the fuel follows 
that given in reference 3. The calculational model reduces the fuel swelling problem to 
a determination of the creep of an internally pressurized thick-walled sphere. The 
multiaxial creep equations of appendix B can be applied to a sphere by making the follow- 
ing substitutions 


°1 " °2 ~ °t 


°3 = a r 


e l ~ e 2 " c t 
e 3 = c r 

Because of the symmetry of components 1 and 2, equations (B5) and (B6) reduce to 

o* = o t - a r 


• * Cx /• • \ 

c = o (e f e r) 


(Cl) 

(C2) 


It is assumed that the volume of the solid does not change during plastic deformation 
(v = 0). This condition of incompressibility reduces equation (B2) to 


Combining equations (C2) and (C3), 


2e t = -^ r 


e* = 2k, 


The strain rates in spherical geometry are given by 


(C3) 


(C4) 


= du 

(C5) 

r dr 


• u 

(C6) 
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where u is the rate at which an incremental radius dr is changing. 

Substituting from equations (C3) and (C6) into equation (C5) and integrating, 



Combining equations (C4), (C6), and (C7), 



The equilibrium equation for a sphere is 

dcr 

r _,2 ( a t -<T r ) 

Combining equations (B4), (Cl), (C8), and (C9), 

"'r . 2 / 2C lV /n r -8/(n-l) 

Integrating equation (CIO), 



The boundary condition at the inside and outside surfaces of the sphere are 


cr . = -d. 
n 


a = -p 
ro F o 


Applying these boundary conditions to equation (Cll), 
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(C7) 


(C8) 


(C9) 


(CIO) 


(Cll) 


111 nnin 


l in lit 


ii n m r 



f 


4 


■t 


» 


♦ 



Differentiating equation (C12), 



Combining equations (B4), (C2), (C4), (C9), and (C13), 



(C12) 


(C13) 


Then at r = r . 



(C14) 


The radii in equation (C14) can be written in terms of volume ratios. Let x^ denote the 
fractional fuel volume change caused by gaseous fission products, that is, the gas vol- 
ume expressed as a fraction of the original fuel volume. Then, 
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X 


(C15) 



and 


dt dt 


(C16) 


The solid fuel volume is assumed here to be constant with respect to time and equal to 
the original fuel volume. Therefore, 


M r °- r 0 =v ° 


(C17) 


and combining with equation (C15) yields 


r. x„ 
_1 g_ 


.3 1 + x„ 


(C18) 


Substituting equations (C17) and (C18) into equation (C14) yields 


1/n 


, 1 + x„ 


n 


dx 


g = 3a 
2 


3 A(i.V 

2n \ 


dt 


(C19) 


Assuming the gaseous fission products to be ideal gases 


Pi = 


2C y KT f Nbt 

X_T 


(C20) 


Combining equations (C19) and (C20) 
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(C21) 



A pressure p Q at the outside surface of the sphere could be produced as a result of 
a constraint at the outside of the fuel body. In order to simplify the problem and yield a 
conservative design, the pressure p Q is taken to be zero. Equation (C21) then becomes 


^3C y KT f Nbt\ n 
nr 


x„ = 
g 


n-!| 
g 



in 


(C22) 


Equation (C22) will therefore yield a volumetric swelling rate corresponding to an 
unrestrained condition which will be greater than the expected swelling rate if a clad is 
present to make p Q greater than zero. 

Integrating equation (C22) between the limits from zero to and zero to t, for a 
constant fission rate b/ r, yields 



Equation (C23) is equivalent to the results of reference 3, except that in reference 3, x 

& 

was assumed to be small compared to 1 and the term (1 + x ) does not appear. 

D 

The swelling caused by solid fission products is assumed, as in reference 2, to be 
directly proportional to the number of fissions. This assumption is based on atomic or 
molecular volumes originally proposed in reference 8. Solid fission product swelling is 
then expressed as 


x 


s 


C s bt 

r 


(C24) 
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where the constant C_ is the fractional change in atomic or molecular volume of solids 

s 

occurring as a result of fission. The value of C depends on the fuel compound. 

S 

The total change in fuel volume expressed as a fraction of the original volume is 

x = x g + x g (C25) 

The engineering bulk strain rate is then 

x = Xg + x g (C26) 

where for a constant fission rate x is given by equation (C22) and x is obtained by 

& s 

differentiating equation (C24) with respect to time. 

* s = C B £ (C27) 

T 

The true bulk strain rate may then be determined from equation (B18). Substituting 
equations (C22), (C25), (C26), and (C27) into equation (B18) yields 


r 


3a /SCKTjNbtV 


nr 


v=< 


, n-1 


1 - 


1/n 


L 1 + x„ 


n 


+ C s-V 


i + x g + 


C bt 


(C28) 


j 


Equation (C28) gives the true bulk strain rate of the fuel for a constant fission rate 
based on the calculational model assumed herein. 


i 


i 
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APPENDIX D 


CREEP ANALYSIS FOR FUEL AND CLAD CYLINDERS 

Stress and strain in the fuel and clad cylinders can be analyzed using the general 
equations of appendix B and the volumetric fuel swelling equations of appendix C. Equa- 
tions (B2), (B5), (B6), and (B8) can be written in polar coordinates for application to 
cylindrical geometry as follows 


e z + + e r “ ^ 

(Dl) 


(D2) 

°* * ^ [(°a - ff t f * (' \ - °r? + (°r ' 

(D3) 


(D4a) 

• 0 *n-l 

^t = 3 + ~^~2 (2a t ” a r ” a z } 

(D4b) 

• *n - 1 

^r = f + i V-( 2 -r- ff z--t) 

(D4c) 


It is assumed that the cylinder is long compared to its diameter so that plane strain 
exists in the axial direction. Therefore, e_ is independent of the radius. Consequently, 
de z /dr = 0 and e z = Cj. 

Applying this restriction, equations (Dl), (D2), and (D4a) become 


C 1 + e t + ^r “ v 


** = y- [( c > ' ^ + & ' ' e 'f + ^ ' c ') 2 ] 


1/2 


(D5) 

(D6) 


27 


Cj = I + 


*n-l 


3 2 


(2a z - a t - <r p ) 


(D7) 


The compatibility equation (ref. 5) is expressed as 


<*t . . 

r = e r " e t 

dr r 1 


(D8) 


Combining equations (D5) and (D8) and integrating, 


e t = 


* -Ci Cg 

2 'r 2 


Let 


<P = — 
r^ 
m 


(D9) 


where r is a reference radius. Then, 
m 7 


v-Ci C 2 


2 2 
V 


(DIO) 


Combining equations (D5) and (DIO), 


e„ = 


v - Ci C 2 
+ • 


2 „ 2 n 2 

W 


(Dll) 


Combining equations (D6), (DIO), and (Dll), 


2C 0 

.* ^ 


* 4 / 


(l + A 2 * 4 )' 


1/2 


(D12) 


where A is defined as 
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(D13) 


A a 


(3C 1 - v)r‘ 

2V3 Co 


Combining equations (B4) and (D12) 


( 7 * = 


'2C 2 Vl + A 2 / 


v l/n 


,/T 2 2„ 

* 3 r m <P a 


(D14) 


The equilibrium equation (ref. 5) is 


dcT r 

r = or* - or 

dr 1 r 


or using the definition of <p (eq. (D9)) 


da 


(D15) 


Combining equations (D4b), (D4c), and (D15), 


da 


cp 


d cp 


3aa 


*n-l 


(e t - e r ) 


Combining this equation with equations (DIO), (Dll), and (D14) gives the radial 
stress gradient 


dor 

d<p 


-B 


(l + A V 4 ) <n ' 1)/2n ?. (n+2)/n 


where B is defined as 


2 n+1 Cr 


,1/n 


B = 


i3 (n + l)/2 a 2 
m> 


(Die) 


(D17) 
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The radial stress can be determined from equation (D16) and the boundary condi- 
tions which are discussed later in this appendix. 

The tangential stress may be obtained from the equilibrium equation (eq. (D15)) and 
equation (D16) 


°t = °r “ 


/ 2 4\( n- l)/2 n 2/n 

(l + A > 4 ) <p Z/n 


(D18) 


The axial stress may be obtained by combining equations (D17), (D14), and (D18) 
using the definitions of A and B (eqs. (D13) and (D17)) 


cr = a + 
z r 


b(aV3 cp 2 - l) 


/ 2 4\( n “l)/2 n 2/n 

2(l + A > 4 ) cp 2/n 


(D19) 


The tangential and radial strain rates are given by equations (Dl) and (Dll) in com- 
bination with equations (D13) and (D17) 


• _ v a3 (n+1)/2 B n / , 

3 2 < n+ V \ Vi 


(D20a) 


i A. 

r 3 2 n+1 <^ 2 \ Vi 


(D20b) 


Boundary Conditions 

The boundary conditions can be specified in terms of the radial stress at the inside 
and outside surfaces of the fuel and clad cylinders, the axial force on the cylinders, and 
the tangential strain rate at the interface of the fuel and clad. 

The pressure at the inside surface of the fuel cylinder and the outside surface of the 
clad cylinder is taken to be zero. Also the fuel and clad cylinders have a common con- 
tact pressure p m at their interface. Therefore, from the definition of cp (eq. (D9)) 

r a 

a = 0 at (p = = <p 

r 
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r, 

a = 0 a t <p = = cp h 

r m 

CT r = "Pm at <P = 1 

where r & is the inner radius of the fuel cylinder, r^ is the outer radius of the clad 
cylinder, and the reference radius r m is taken to be the radius at the fuel-clad inter- 
face. Applying these conditions to equation (D16) and integrating yields 


= B. 


a 


r 




/ n 

n b 

J(p 

f 

J( Pu 


dcp 


(X ♦ &f- 1)/2 y»+»/n 


d <p 


(l ♦ A^«) <n ' l>/2 V n+2>/n 


dcp 


(l ♦ Af V) (n - 1,/2 V-’)/n 




&<P 

(i ♦ a 


(D21) 


(D22) 


(D23) 


(D24) 


It is assumed that no axial force is transmitted between the fuel and clad cylinders. 
Therefore, the net axial force on the cylinders is zero, or 



a z cp dcp = 0 



a z cp dcp = 0 
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Applying these conditions to equations (D19), (D21), and (D22), 


F f F* 

jy j<p 


dip 


A f Vi cp 2 - 1 


(i + Ayf' 1)/2n ^ 2 ^ 2(1 + a f v ) (n ' 1)/2 V /n 


cp dcp-0 (D25) 


/[/ 




d<p 


A V3 cp 2 - 1 


(i + Ay) <n ' 1)/2 V n+2)/n 2(1 ♦ A y) (n - 1>/2 V/" 


(p d cp = 0 (D26) 


The tangential strain rates of the fuel and clad must be equal at their interface be- 
cause they are in contact. Therefore, applying equation (D20) to the fuel and clad cyl- 
inders at their interface (cp = 1), 




a 3 
c 


(n c +l)/2 


~tm 


n +1 
1 c 



(D27) 


(D28) 


When applying equation (D20) to the clad, v was taken to be zero. Also, it must be 
recognized that equation (B4), which was utilized in the derivation of equation (D20), is 
an empirical equation and does not specify the absolute sign of the strain rate. There- 
fore, the stress must be inserted as a positive value and the strain rate specified as 
positive or negative depending on whether the stress is positive (tensile) or negative 
(compressive). 

Combining equations (D27) and (D28) to eliminate the pressure p m at the interface 
yields 


t 
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3e 


tm 


= 1 - 


V n c 


^.vy-v 


v\ 1 - ( V n c> 


© 


^x(Vn f )-(V4 , c ) 



in. 


(D29) 


or 


where 


and 



(D30) 


(D31) 


(D32a) 


(D32b) 


The quantity /3 is a function only of n and the radius ratio <p a or <p b> The latter is 
evident from equations (D23) and (D25) or (D24) and (D26). Values of j3 were computed 
by numerical integration of equations (D25) and (D26) to determine A and by numerical 
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integration of equations (D23) and (D24) to determine P m /B. The results of these com- 
putations for j3 and A are shown in figures 2 and 3. 

The tangential strain rate at the fuel-clad interface can be obtained from equations 
(D30) and (D31), and figure 2. All other strain rates can be written in terms of the tan- 
gential strain rate at the interface. 

Equation (D20a) may be rewritten as 



(D33) 


Equation (D20b) may be rewritten as 



Then from equation (D5) 



(D34) 


(D35) 


When applying equations (D33), (D34), and (D35) to the clad v is taken to be zero. 
The strain rates are a function of time t and must be integrated with respect to 
time to obtain the strain in the time period from time zero to time t as follows 


e = 



(D36) 


Equations for radial, tangential, and axial stress can be rewritten in terms of the 
constants evaluated from the boundary conditions. Solving for the radial stress at the 
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fuel-clad interface p m from equation (D28) and substituting for B/p m from equation 
(D32b) 



(D37) 


The tangential and axial stresses at the inner and outer surfaces of the fuel and clad 
cylinders relative to p m are obtained from equations (D18) and (D19) using the definition 
of 0 (eqs. (D32a) and (D32b)): 



L a 2 4\ (n f -1) / 2n f 2/ n f[Vi L 
^1 + Af (p a j <P a — (1 
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TABLE I. - SELECTED DATA AND RESULTS FOR NUMERICAL EXAMPLE 


Time, 

t, 

hr 

Ratio of 
time to 
total 
irradia- 
tion 
time to 
achieve 
fraction 
b, 
t/r 

Fractional 
fuel volume 
change 
caused by 
gaseous 
fission 
products, 

x g 

True bulk 
strain rate 
caused by 
solid and 
gaseous 
fission 
products, 

V 

Quantity 

defined 

by 

equation 

( D 31 ), 

a 

Relative 
tangen- 
tial 
strain 
rate at 
fuel- 
clad 
inter- 
face, 

3 «tm^ 

Tangential 
strain rate 
at fuel- clad 
interface, 

e tm 

Tangen- 
tial 
strain 
at fuel- 
clad 
inter- 
face, 

6 tm 

| 

Tangen- 
tial 
strain 
at inner 
radius 
of fuel 
cylinder, 

6 ta 

Tangen- 
tial 
strain 
at outer 
radius 
of clad 
cylinder, 

6 tb 

Axial 

fuel 

strain, 

6 zf 

Axial 

clad 

strain, 

^zc 

0 

0 

0 

0 . 180 X 10" 5 

0.628 

0.626 

0 . 376 X 10" 6 

1 

0 

0 

0 

0 

0 

5 000 

.25 

.021 

.761 

.708 

.599 

1.52 

.006 

-.004 

.004 

.016 

-.002 

10 000 

.50 

.059 

1.02 

.727 

.593 

2.02 

.015 

-.010 

.011 

.041 

-.006 

15 000 

.75 

.114 

1.28 

.740 

.589 

2.50 

.026 

-.019 

.020 

.073 

-.010 

, 20 000 

1.00 

oo 

CO 

1.55 

.752 

.585 

j 

, 3.01 

.041 

-.030 

.030 

.114 

-.015 


r c\i 

/ 

'/mm 

7M//////// 

,i/g 

wmm 


Figure 1. - Schematic drawing of cylindrical fuel element (fuel pin) showing 
fuel, clad, and void space. 
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(b) Clad cylinder. 

Figure 2. -The quantity $ as function of 
radius ratio and n for fuel and clad cylin- 
ders, for use in equation (D31). 


* n < 

tnmmi 



-ii 



TO 

.i tt- "rjT 

-ii 




“ l;il -iii 


_ / il p 


!■• ; jj lit 

'll! . l " 
i't. ■ . j ! iiL 

ji ; . i j i 


Virt 

!|.i iv ! , 

|jj j p 

, |; • !' ]r 

':!! ii’ iir 



r 


(a) Fuel cyclinder. 


* 'l 

,ii iii! 1 ii ' V ' S.'« '•••K 1 ? 

il. iii'-:' ;:•':■!! : ,!i. il * i.jil 

iii- : II!' »! ' ' 'ii . !i! i 


*1.0 1.4 1.8 2.2 2.6 3. 

Radius ratio, <p b 

(b) Clad cylinder. 

Figure 3. - The quantity A as function of radius ratio and 
n for use in equations (D33), (D34), and (D35). 
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Figure 4. - Graphical presentation of equation <C23). Volumet- 
ric fuel swelling caused by gaseous fission products as func- 
tion of parameter defining fuel characteristics, design condi- 
tions, and operating time. 



Figure 5. - Graphical presentation of equation (D30). Relative tangential strain rate at fuel-clad 
interface as function of parameter a defined by equation (D31). 
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Engineering strain, 6 True strain rate. « tm , hr' 



Operating time, t, hr 


Figure 6. -True tangential strain rate at fuel-clad interface as function of operating time 
for numerical example. 



Operating time, t, hr 

Figure 7. - Tangential and axial strains in fuel and clad as function of operating time for 
numerical example. 
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